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Abstract. Using the negative binomial distribution (NBD) and the generalized Glauber-Lachs (GGL) 
formula, we analyze the data on charged multiplicity distributions with pseudo-rapidity cutoffs ij c at 0.9, 
2.36, and 7 TeV by ALICE Collaboration and at 0.2, 0.54, and 0.9 TeV by UA5 Collaboration. We confirm 
that the KNO scaling holds among the multiplicity distributions with r) c = 0.5 at yfs = 0.2~2.36 TeV and 
estimate the energy dependence of a parameter 1/k in NBD and parameters 1/k and 7 (the ratio of the 
average value of the coherent hadrons to that of the chaotic hadrons) in the GGL formula. Using empirical 
formulae for the parameters 1/k and 7 in the GGL formula, we predict the multiplicity distributions with 
r/ c = 0.5 at 7 and 14 TeV. Data on the 2nd order Bose-Einstein correlations (BEC) at 0.9 TeV by ALICE 
Collaboration and 0.9 and 2.36 TeV by CMS Collaboration are also analyzed based on the GGL formula. 
Prediction for the 3rd order BEC at 0.9 and 2.36 TeV are presented. Moreover, the information entropy is 
discussed. 

PACS. 13.85.Hd - 



1 Introduction 



Very recently ALICE Collaboration has investigated the 
multiplicity distributions with pseudo-rapidity cutoffs (rj c = 
0.5, 1.0, and 1.3) at ^/i = 0.9 and 2.36 TeV p. Therein, 
it has compared its data with the data at 0.2, 0.54, and 
0.9 TeV by UA5 Collaboration [1[3], and concluded that 
the combined data with r\ c — 0.5 at 0.2, 0.9, and 2.35 
TeV are fairly well described by the single NBD (negative 
binomial distribution) [H[S]. Moreover, ALICE Collabora- 
tion has reported that the KNO scaling [6j holds among 
the combined data with r\ c — 0.5 at 0.2, 0.9, and 2.35 
TeV. The first aim of this paper is to confirm the state- 
ment above mentioned in pQ and to analyze the same data 
by the GGL (generalized Glauber-Lachs) formula [7J?]. To 
investigate possibility of predictions on multiplicity dis- 
tributions with rj c = 0.5 at 7 and 14 TeV [9 , empirical 
formulae (for 1/k and 7) are adopted. 

Moreover, ALICE and CMS Collaborations have re- 
ported the data on Bose-Einstein correlations (BEC) at 
0.9 and 2.36 TeV [TOlfTTj . Thus we are going to inves- 
tigated them based on a conventional formula with the 
degree of coherence and the GGL formula. 
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First of all, the NBD is introduced in the following 

r(n + k) ((n)/k) n 



P k {n) = 



r(n+ l)r(k) (1 + (n)/k) n + k ' 



(1) 



where (n) and k are the average multiplicity and the in- 
trinsic parameter, respectively. In the KNO scaling limit 
(n and (n) are large, but the ratio z = n/(n) is finite), for 
the quantity (n)P(n, (n)) the following gamma distribu- 
tion is derived from Eq. ([T]) as 



r(k) 



z k-l e -kz 



(2) 



Second we turn to the GGL formula which is expressed as 
follows: 



Pk(n) = 



(p(n)/k) n 



(1 +p(n}/k) n+k 
. L {k-l)f 



exp 



iv{n) 



1 + p(n) jk 



1 + p(n) jk 



(3) 



where 7 = |C| 2 /^4 (the ratio of the ratio of the aver- 
age value of the coherent hadrons to that of the chaotic 
hadrons), p = 1/(1 + 7), and L^ 1 ^ stands for the La- 
guerre polynomials, respectively. Here it should be stressed 
that the GGL formula has three limits, the original GL 
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formula (k = 1) [H], the NBD (7 = 0), and the Poisson 
distribution (fc = 1 and 7 = 00, or 7 — and fc = 00). 

(k = 1) original GL 
Eq. © - (7 = 0) NBD 

(fc = 00) 

(fc = 1, 7 = 00) Poisson distribution 

The KNO scaling function of Eq. ([3]) is given in the 
following 



ipk(z, p) 



■ exp 



(1-p + z) 

P 



4-i 



fc-i 



2y/z(k/p)*(l-p))(4) 



where Ik-i is the modified Bessel function. Eq. © be- 
comes the gamma distribution, as 7 = 0. 

In order to analyze of Bose-Einstein correlations (GGLP 
effect [13], or hadronic HBT effect [HQ23) at LHC, we are 
going to use the following formulae: The first one is well 
known as the conventional formula, 

AT( — ) /jV BG (conventional formula) = 1 + A B| s , (5) 

where the parameter A is named the degree of coherence, 
and E2B is function of momentum transfer (Q 2 = —(pi — 
P2) 2 ) and the range of interaction R. E 2 b = exp(— R 2 Q 2 ) 
(Gaussian formula) and/or E 2 b — exp{—R\/Q 2 ~) (expo- 
nential formula) are used. The second one proposed in 
Ref. [14] is relating to the GGL formula as follows 



and 2.36 TeV. Hereafter we use the CERN-MINUIT pro- 
gram. Results are shown in Fig. Q] and Tabic [T] From re- 
sults based on analyses of the NBD, we observe that 1/fc 
increases gradually as y/s increases. On the other hand, 
the estimated sets of (1/fc and 7 ) in the GGL formula 
show different behavior. Using the results of parameters 
for data with r\ c — 0.5, in particular, we estimate the en- 
ergy dependence of 1/fc in the NBD, and those of 1/fc and 
7 in the GGL formula as follows: 

1/fc (NBD) = Q49 + . 12 ln(V5/0.2), (7) 
!/ fc (GGL) = g9 + Q31 ln (^/ .2), ( 8 ) 
7 (md) = 2 1Q _ .551n(yi/0.2), (9) 

where MD stands for the multiplicity distribution. Energy 
dependences of parameters l/fc (NBD) , l/fc (GGL) and 7^°) 
are shown in Fig. [5J 

Moreover, the data with r\ c = 1.0 at 7 TeV have been 
reported by ALICE Collaboration [5] . The data are relat- 
ing to inelastic events with a positive condition INEL> 
and no value on (n), nevertheless. Then in our analyses, 
the following modifications with a normalization factor c 
is adopted: 



>/iV SG ((GGL)) = 1 

1 + [2p(l 



[2p{l - p)E 2B 

P)E 2 B 



P 2 E 2 2E 



\/k 

P 2 E 2 2B ], (6) 



where p = 1/(1 + 7), and (fc — » 1) means the identical 
charged ensemble. 

The present paper is organized in the following: In 
the second paragraph, we analyze the data with pseudo- 
rapidity cutoffs by means of the NBD and the GGL for- 
mula. In the third paragraph, the distributions of the KNO 
scaling are analyzed by Eqs. @ and (Ql. In the fourth 
paragraph, we consider the physical meaning of the pa- 
rameter 7 in the GGL formula: Data on the 2nd order 
BEC at 0.9 and 2.36 TeV [TU1ITT] are analyzed by Eqs. © 
and ([6]) 14. 15 . In last paragraph, concluding remarks and 
discussion for the information entropy are given. 



2 Analyses of data on multiplicity 
distributions by the NBD and the GGL 
formula 

Utilizing Eqs. ((T|) and ([3]), we analyze the data with pseudo- 
rapidity cutoffs (r/ c = 0.5, 1.0, and 1.3) at 0.2, 0.54, 0.9 



Eq. (P) -> c x Eq. (Q} (with free (n)), 
Eq. @ ^ c x Eq. © (with free (n)), 



(10) 
(11) 



The estimated values are added in Table \T\ It can be said 
that multiplicity distribution at 7 TeV is fairly well ex- 
plained by the GGL formula as well as NBD. (See Fig.[TJ) 
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Fig. 2. Energy dependence of parameters l/fc (NBD) , l/fc (GGL) , 
and 7 (MD) for data with Vc = 0.5. Values at 0.9 TeV by UA5 
Collaboration are omitted in GGL formula, because of extreme 
error bars. 



To get more useful physical information on multiplic- 
ity distributions, the KNO scaling distributions are inves- 
tigated in the next paragraph. 
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Table 1. Results in analyses of multiplicity distributions (NSD except for 7 TeV). 
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UA5, = 


0.2 TeV 






0.5 


0.56±0.03 


11.5/13 


1.00±0.22 


2.2±0.2 


3.8/12 


1.0 


a a cr_i_n no 
0.4o±0.02 


i n /oo 
19/32 


n on_i_n Ad 
O.o9±0.4o 


2.0±2.2 


10. ( /31 






UA5, ys = 


U.54 lev 






0.5 


0.59±0.02 


30/22 


0.86±0.19 


1.4±0.7 


21/21 


1.5 


0.50±0.01 


13.8/28 


0.52±0.03 


0.22±0.26 


13.5/27 






UA5, vs = 


0.9 TeV 






0.5 


0.65±0.03 


2.1/22 


0.65±0.03 


0.0±11.3 


2.1/21 


1.5 


0.56±0.02 


18.7/51 


0.56±0.02 


0.0±63.0 


18.7/50 






ALICE, 


= 0.9 TeV 






0.5 


0.69±0.01 


10.4/23 


0.89±0.13 


0.88±0.42 


3.7/22 


1.0 


0.65±0.01 


44/41 


0.79±0.08 


0.74±0.28 


35/40 


1.3 


0.61±0.01 


70/47 


0.74±0.08 


0.70±0.29 


62/46 






ALICE, 75 = 


= 2.36 TeV 






0.5 


0.82±0.02 


19/26 


l.OOiO.ll 


0.75±0.08 


14/25 


1.0 


0.73±0.02 


82/45 


0.79±0.09 


0.41±0.33 


80/44 


1.3 


0.67±0.02 


145/51 


0.70±0.04 


0.25±0.20 


144/50 






ALICE, y/s = 7 TeV (inelastic, n 


>0) 




1.0 


0.88±0.01 


217/62 


0.99±0.03 


0.45±0.08 


191/61 




(n) = 11.6±0.1, c= 1.07±0.01 


(n) = 11.5±0.1, c= 1.08±0.01 



Table 3. Results in analyses of KNO scaling distributions (n)P(n) described by KNO variable z = n/{n). 





NBD 




GGL 




Vc 


l/k 


X 2 /NDF 


l/k 


7 


X 2 /NDF 






UA5, . 


/s = 0.2 TeV 






0.5 


0.92±0.04 


59/12 


l.OOiO.01 


0.72±0.12 


38/11 


1.5 


0.53±0.02 


28/31 


1.00±0.13 


2.2±0.1 


14/30 






UA5, x 


fs = 0.54 TeV 






0.5 


0.87±0.02 


172/21 


l.OOiO.OO 


0.77±0.06 


100/20 


1.5 


0.58±0.01 


41/28 


0.75±0.04 


0.93±0.18 


24/27 






UA5, - 


/i = 0.9 TeV 






0.5 


0.91±0.04 


22/21 


1.00±0.02 


0.59±0.11 


13/20 


1.5 


0.62±0.02 


15.3/50 


0.62±0.02 


0.0±17.0 


15.3/49 






ALICE, 


= 0.9 TeV 




0.5 


0.92±0.01 


80/22 


l.OOiO.01 


0.44±0.04 


54/21 


1.0 


0.76±0.01 


99/40 


l.OOiO.01 


0.93±0.04 


56/39 


1.3 


0.69±0.01 


110/46 


1.00±0.03 


1.20±0.04 


73/45 






ALICE, 


= 2.36 TeV 




0.5 


1.00±0.05 


94/25 


l.OOiO.01 


0.31±0.08 


87/24 


1.0 


0.79±0.02 


128/44 


1.00±0.03 


0.85±0.06 


118/43 


1.3 


0.73±0.02 


159/50 


1.00±0.09 


1.07±0.06 


156/49 


combined UA5, 0.2, 0.9 TeV, ALICE, 0.9, 2.36 TeV 


0.5 


0.94±0.01 


264/83 


l.OOiO.OO 


0.44±0.03 


203/82 



3 Analyses of data on KNO scaling 
distributions by Eqs. (J2D and (Sj) 

Utilizing the KNO scaling variable z(= n/(n)), data on 
the KNO scaling distributions (n)P(n, (n)) are shown in 
Fig. [3J Our results analyzed by Eqs. ^ and (0} are given 
in Table H 

Second, we combine the data with rj c = 0.5 at 0.2, 0.9 
and 2.36 TeV and analyze them by Eq. © (the gamma 



distribution) and Eq. Q (the modified Bessel function). 
Since data with r\ c = 0.5 at 0.54 TeV are not well explained 
by two formulae, three data are regarded as exceptional 
among them. In the final paragraph, we consider this fact. 
Results on data with ry c = 1.3 at 2.36 TeV show large x 2 
values; Large values of x 2 's (x 2 >> NDF) denote that the 
single NBD or the single GGL formula should be improved 
for explanations of data. See, for example, Ref. |16) . 
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Table 2. Energy dependence of parameters (l/fc' NBD \ 
l/fc (GGL) , and 7 (MD) ) for multiplicity distributions with t) c = 
0.5. Eqs. ©, ©, and © are used. See Fig. [2] 





l/k^ NBD > l/k^ L > 


7 (MD) 


0.9 


0.69 0.89 


0.88 


2.36 


0.82 1.00 


0.75 




calculated values 




7 


0.93 1.0 


0.14 


14 


1.00 1.0 


0.0 



GeV/c. Using estimated values of <y( BEC ) and R in the 
2nd BEC by ALICE Collaboration, we can predict the 
3rd order BEC; The following formula [2] is used, 

N {z-)/ N bg = i + 6 p(l-p)exp (-lRyjQ%\ 



(15) 



+3^(3-2p)exp^--i?\/Q§ 
+2p 3 cxp (-RJC%) , 



Moreover, the results on 1/k in Table [3] by Eq. ^ 
(gamma distribution) show similar behavior to those in 
Table [TJ On the other hand, the sets of (1/k and 7) esti- 
mated by Eq. (U|) are almost approximately (1/k = 1.0) 
at 0.2~2.36 TeV. The parameter 7 depends on colliding 
energies, which is roughly expressed as 



where p = 1/(1 + 7) an d Q§ = Q12 



7 



(KNO) 



0.96 - 0.29^(^/0.2) 



(12) 



Comparing Eq. (fT2"|) with ©, we observe that 7 



(KNO) 



7 



(MD) 



/3. From the bottom panel of Fig. [3j we know that 
coincidence among combined data with r\ c = 0.5 (the num- 
ber of data points is 84) is fairly well explained by the GGL 
formula with smaller value of % 2 than that of the NBD. 



4 Analyses of data on the 2nd order BEC by 
means of Eqs. ([5]) and ([5]) 

We analyze the data on BEC at LHC by the use of Eqs. (|5|) 

and (© with E 2 b — exp(— R 2 Q 2 ) and/or E 2 b = exp(— Ry 
The following modification, i.e., introducing the normal- 
ization factor, is used, c is reflecting the long range corre- 
lation (1 + SQ)" 1 in many data. 

Q- 



Eq. © 
Eq. © 



c x Eq. ©- 
c x Eq. ©- 



->c(l 



1 



A), 
(1 



2 7 



(BEC) 



(1 + 7 (BEC))2 



(13) 

(14) 



Results are depicted in Table |4] and Fig.|U In Eq. (0, tne 
effective degree of coherence "A" is "(1 + 27)/(l + 7) 2 ". 
In our concrete analyses, we obtained that l/fc' BEC ) = 1. 
Because estimated value of this seems to be reasonable 
for the identical charged ensemble, it is omitted. 7( BEC ) is 
similar to the value at 0.9 TeV by ALICE Collaboration 
in Tabled! 

Furthermore, by the use of Eqs. © and (O with c, we 
have analyzed the data on BEC at 0.9 and 2.36 TeV by 
CMS Collaboration [TT]. Results are shown in Fig. @] and 
Tabled) R is emphasized that the ratio 7( BEC ) decreases, 
as the colliding energy increases. In other words, the effec- 
tive degree of coherence "A" and the range of interaction 
R increases from 0.9 to 2.36 TeV. To draw more signifi- 
cant meaning about the parameter 7, we need the BEC 
measured in data with r] c — 0.5. 

It is worthwhile to predict the 3rd BEC at 0.9 TeV 
using the same condition with M < 6, 0.1 < fcx < 0.55 



Qh + Qii- Our 

prediction is depicted in Fig. @|d). 

Furthermore, predictions on the 3rd order BEC at 0.9 
and 2.36 TeV for CMS Collaboration are also displayed in 
Fig-Hd) and (e). 
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60 



Fig. 5. Expected multiplicity distributions with r/ c — 0.5 
at t/s — 7 and 14 TeV. Computations are based on the 
GGL formula (Eq. ©) with values in Table and (n) = 
2.5 + 0.76 ln(vV0.2). 



5 Concluding remarks and discussion 

Through our present analyses, we summarize our conclud- 
ing remarks (C1~C4), and add discussion (Dl) for the 
information entropy. 

CI We have confirmed that the multiplicity distributions 
with T] c = 0.5 are described by the single NBD pQ. 
Moreover, we also confirm that the GGL formula does 
work well for the explanation of the same data in present 
analyses. As the pseudo-rapidity cutoffs increase, the 
i] c = 1.0, and 1.3 at y/s = 0.2, 0.9 and 2.36 TeV show 
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Table 4. Analysis of data on BEC by ALICE Collaboration and CMS Collaboration. Because estimated value of l/fc' BEC ' is 
a unit, it is not cited. 



Eq. H3J | Eq. JUJ 

(upper: Gaussian formula, and lower: exponential formula) 
X c R (fm) x 2 /NDF | 7 c R (fm) X 2 /NDF 



■s/s = 0.9 TeV, ALICE (multiplicity M < 6, 0.1 < k T < 0.55 GeV) 


0.35±0.02 0.99±0.00 0.83±0.04 121/72 
0.64±0.04 0.98±0.00 1.33±0.09 98/72 


4.0± 0.3 0.99±0.00 0.81±0.03 119/72 
1.30±0.23 0.98±0.00 1.18±0.07 98/72 


v/s = 0.9 TeV, CMS (Excluding 0.6 < Q < 0.9 GeV/c) 


0.32±0.01 0.99±0.00 0.96±0.02 407/165 
0.66±0.02 0.99±0.00 1.75±0.04 229/165 


4.5± 0.2 0.99±0.00 0.95±0.02 394/165 
1.08±0.13 0.99±0.00 1.59±0.03 225/165 


v/s = 2.36 TeV, CMS (Excluding 0.6 < Q < 0.9 GcV/c) 


0.33±0.03 1.00±0.00 1.20±0.07 80/81 
0.72±0.08 1.00±0.00 2.32±0.17 75/81 


4.3± 0.6 l.OOiO.OO 1.18±0.07 80/81 
0.84±0.39 1.00±0.00 2.03±0.08 76/81 



slightly weak violations in KNO scaling distributions, 
because of large values of (x 2 /NDF)'s. 

C2 Estimated x 2 m the GGL formula are slightly better 
than those of the NBD. In Table [3J we observed that 
distributions with r\ c = 0.5 at 7 TeV does not have the 
coherent component. In other words, the multiplicity 
distributions with r\ c — 0.5 at 7 TeV are described by 
the NBD with k = 1. 

C3 Using values in Table [21 we can predict multiplicity 
distributions with r\ c — 0.5 at 7 and 14 TeV in Fig. [SJ 
Those are able to be examined in a near future. If 
there were discrepancies among data and predictions, 
we should consider the other effect, for example, due 
to the mini-jets |16j . 

C4 Through present analyses of the BEC, results by the 
exponential formula seem to be better than those by 
the Gaussian formula in Table|H See [T7] for the source 
functions. Moreover, values of 7's obtained in Tables Q] 
and |4] seem to be similar each other. To obtain more 
significant knowledge on the parameter 7, analyses of 
the multiplicity distributions and the BEC in the same 
hadronic ensembles are necessarv |14U15) . Our predic- 
tions for the 3rd order BEC at s/s = 0.9 and 2.36 GeV 
would be compared by measurements as UA1 Mini- 
mum Bias Collaboration did [18] . By the comparisons, 
we could obtain more useful information on the pa- 
rameter 7 and the role of the GGL formula. 

Dl We have to consider reasons of the large value of x 2 
concerning data with rj c = 0.5 at 0.54 TeV in Table |3l 
For this aim, we use the information entropy defined 
as [TO], 



S , = -^P(n)lnP(n). (16) 



Results are depicted in Table [5] From them, we calcu- 
late a new plot of (t] c /Y max , S/Y max ), where Y max = 
hi^/s/mp. We see a kind of scaling law on the infor- 
mation entropy for multiplicity distribution [3]. 



Table 5. The information entropy of data by ALICE Collab- 
oration and UA5 Collaboration. AS — — ^2 8P(n) In P(n) — 
Y J 5P{n). Theoretical values are S (NW) S (GGL) ^ S± 0.03. 



(TeV) 


7? c S AS 


0.2 


0.5 2.05 0.10 
1.5 3.03 0.25 


0.54 


0.5 2.23 0.07 
1.5 3.22 0.13 


0.9 (UA5) 


0.5 2.37 0.13 
1.5 3.37 0.29 


0.9 (ALICE) 


0.5 2.39 0.11 
1.0 3.03 0.18 
1.3 3.27 0.22 


2.36 


0.5 2.56 0.15 
1.0 3.22 0.25 
1.3 3.41 0.29 


7 


1.0 3.45 0.08 



Next, we consider the information entropy in the KNO 
scaling, which is calculated by the following formula, 

5 (KNO) - ln(n) = - J ^j(z) \ni>(z)dz, (17) 

where z — n/(n). Our results at 0.54 and 0.9 TeV 
are given in Table El where theoretical S"s are given 
as s(Eq.©) and s .(Eq.SJ)_ We gee that the differences 

SS = S< data ) - S {thcoiy) between data with rj c = 0.5 
at 0.54 and 0.9 GeV and theoretical values are larger 
than those of other cases. The large values of the ratio 
X 2 /NDF at 0.54 and 0.9 TeV in Table[Hon KNO scal- 
ing are also observed as the large values of t55"s with 
r) c = 0.5. 
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Fig. 1. Analyses of data with \r]\ < r\ a by means of Eqs. ([T]) and ©. 
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Fig. 4. Analysis of data on BEC at 0.9 TeV by ALICE Collaboration with conditions M < 6, and 0.1 < k T < 0.55 GeV and at 
0.9 and 2.36 TeV by CMS Collaboration. Panels (d) and (e) are also our predictions of the 3rd order BEC for at 0.9 and 2.36 
Te. 



